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^ ■ Abstract. The present work is a generalisation of the dynamical screening factor 

^ , presented in [1] to consider an atom located at an arbitrary position within the 

fuUerene. A more elaborated investigation into the case where the atom is located 
at the centre is performed and compared with quantum mechanical calculations for 
I dynamical screening factor of ArQCeo [2] and MgdCgo [3]. The tt and a plasmons of 

C/3 ' the fuUerene are accounted for in a modified screening factor to improve correspondence 

with the quantum calculations. The spatial dependence of the screening factor was 
explored with ArQCeo and Ar@C24o and found to depend significantly on the radial 
(-H ■ distance of the atom from the centre of the fuUerene. A spatial averaging of the 

CL|! screening factor is presented. 

^ ! PACS numbers: 32.80.Fb, 36.40.Gk, 36.40.Vz 

1. Introduction 

Particles may be confined by a wide range of metliods. Tfiese include trapping by 
00 ' 

Q I optical means, confinement within materials, e.g. as dopants in solids, or within hollow 
structures such as fullerenes and nanotubes. As the ability to insert atoms or molecules 
into these carbon cages increased, so too has the interest in these endohedral systems 



X 



^ I grown [4, 5] . The possible applications of this novel form of matter ranges from forming 
the register of a quantum computer [6-8] to enhancing imaging techniques [9]. 

Understanding the properties of such systems is therefore of much interest. One of 
the means of achieving this is by spectroscopic means. If the confining effect of the cage 
on the optical properties of the endohedral species can be understood, this may provide 
a means of predicting the properties of the whole system. 

The fuUerene cage screens the endohedral atom from external electromagnetic fields. 
Recent theoretical studies [1,10,11], using semi-classical methods, show that this effect 
is strongly dependent on the photon energy. It was found that in the vicinity of the 
plasmon resonances of the fullerene, the photoabsorption of the confined atom is strongly 
enhanced from that of the free atom. A dynamical screening factor may be defined to 
relate the two photoabsorption cross sections. 
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There are other considerations on how the photoabsorption or photoionization 
of the atom is modified by the confining fuUerene cage. Confinement resonances, as 
discussed in [12, 13], consider the interference caused by the photoelectron refiected 
by the fullerene's inner and outer surfaces. The modification of the cross section 
owing to the interference between the confinement resonances and intra-doublet atomic 
resonances are investigated in [14], with reference to the experimental work of [15]. This 
was elaborated upon in [16, 17]. Reference [18] discusses the effect of multi-electron 
correlation on the confinement resonances. Also, the role of multi-walled fuUerene cages 
is studied in [19]. 

In our previous work on the dynamical screening by the fuUerene [1, 10], the 
endohedral atom is located at the centre of the fuUerene. However, interactions with 
the fullerene, thermal vibrations, electron transfer, etc, leads to the atom being not 
necessarily fixed at any one position and that the dynamical screening factor should be 
considered over all possible positions of the confined species. In [1], the dynamical 
screening factor was presented for the limiting case where the mutual interaction 
between the polarized atom and the polarized fullerene cage was neglected. Within this 
approximation, the screening factor is independent of the position of the endohedral 
atom. In the present paper, a more detailed analysis of the case of the centrally 
positioned atom is performed. Also, a comparison of this model with results from time 
dependant local density approximation calculations on ArtQlCeo [2] and MgQCeo [3] is 
elaborated upon. The tt and a plasmons of the fullerene are accounted for in a modified 
screening factor to improve correspondence with the quantum calculations and, on this 
basis, to achieve a better understanding of the physics of the dynamical screening effect. 
The spatial dependence of the dynamical screening factor is investigated for two case 
studies: ArQCgo and Ar@C24o- A method for spatial averaging and its results are 
presented. 

The atomic system of units, e = m = h = 1 is used throughout the paper. 



2. Theoretical Framework 



When the endohedral system is exposed to an external electromagnetic field, the 
fullerene, being dynamically polarized, will screen the confined atom. Depending on 
the frequency of the light, the field at the atom may be stronger or weaker. Therefore, 
the photoabsorption cross section of the confined atom can differ from that of the free 
atom in the same external field. A dynamical screening factor JF = J^{uj) may be defined 
to relate the two: 

^=^. (1) 

A semi-classical approach is taken here to study this effect. The fullerene is 
modelled as a dielectric shell [20-24] and classical electrostatic methods are used. As 
in [1,11], the fullerene is treated as a spherical shell of some thickness Ai? with an inner 
and outer radius of i?i and i?2 respectively. The dielectric function of the fullerene is 




Figure 1. The position of the point-hke atom, indicated by the filled circle, is confined 
at p inside the fuUerene. The angle between p and the external field Eq is denoted 9. 
The field at the position of the fuUerene is given by Etot (p, ^) ■ 



e = e(^)- Given the strong delocalisation of valence electrons of the fuUerene, they 
are treated as a free electron gas. The dielectric function of the fuUerene e(u;) may be 
expressed via the Drude model [25]: 



The plasma frequency ujp is: 



where N is the number of valence electrons and V = 47r(i?2 ~ -Ri)/3 is the volume of 
the fuUerene shell. 

The point-like endohedral atom is characterised by the dynamic dipole polarizability 
a{ui) and can be located at an arbitrary position p within the fuUerene, see figure [1] 

We consider the case when the endohedral system is in a uniform and monochro- 
matic external field, whose strength is given by Ecxt(^) = Eqc"^*, where Eq is a real 
amplitude. 

The photoabsorption cross section of an object is proportional to its 
photoabsorption rate, Q, which can be calculated from the relation [23,26]: 



Here E(r) and dD(r) are the electric field and the elementary dipole moment at the 
position r. The integration is carried out over the volume of the considered system. For 
a free point-like atom the photoabsorption rate becomes 




(2) 




(4) 



Qfree(cc^) = - Im (E*D) = - E2 Im {a{u)) , 



(5) 
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with D = D(ci;) = a(co')Eo being the induced dipole moment of the atom. Similarly, the 
photoabsorption rate of the confined atom is 

Qconf (p, o;) = I Im (E*„t(p, w)Dtot(p, ^)) = | ^lotiP^ ^) Im , (6) 

where Etot(p, t^) is the total electric field at the position p and Dtot(p, t^) is the atomic 
dipole moment, which may depend on the position of the atom and the frequency of the 
field. Thus, the dynamical screening factor becomes 

-^(P> ^) = 7^ — = ^2 — • (") 

^;free -^o 

In order to calculate the dynamical screening factor, it is necessary to calculate 
the total field at the position of the atom. The total field deviates from the external 
field: Etot = Eg + AE(p, tu). Here, the additional term AE(p, cj) appears as a result 
of the fuUerene's polarization. The latter is due to the action of both the external field 
and the field of the dipole induced at the atom. If the atomic polarizability a{uj) is 
small, the inverse effect of the polarized atom on the shell can be ignored. In this case, 
the term AE(p, uj) is solely due to the external field so that it does not depend on 
the position vector p. This means that the total field and, consequently, the dynamic 
screening factor are not sensitive to the atom's position inside the fuUerene cage. Such 
an approach was adopted in previous studies [1]. 

However, the polarizability of the endohedral atom may be large at certain photon 
energies, primarily in the vicinities of giant resonances in the atomic photoabsorption 
cross section. In these energy regions the inverse effect of the polarized atom on the cage 
may become significant. As a result, the induced moments of the atom and the cage 
become correlated and acquire the dependence on p. To calculate the additional field 
AE(p, cij) one has to carry out a self-consistent consideration which takes into account 
the mutual polarization of the atom and the fullerene. 

The self-consistent method is based on the iterative procedure developed in [1]. For 
the sake of clarity, let us outline the main points of this procedure. 

To begin with, one makes use of the following two auxiliary problems. 

The well-known first problem (see, e.g. [25]) concerns the field in the interior of a 
dielectric spherical shell exposed to the uniform external field Eq. The field in question 
can be written as (1 — 2;)Eo, where the quantity z depends on the dielectric function e 
of the shell, and on its inner and outer radii. In the case of a monochromatic field z 
also acquires the dependence on oj. 

The second problem is the determination of the electric field E due to the 
polarization of the shell by a point-like dipole D located at the position p inside the shell. 
Using standard methods of electrostatics, one can demonstrate that at the position p 
this field can be represented as follows: 

E(p) = 5(p)D = --1 [.i(p)D + s^{p) (Dn) n] . (8) 

Here n is the unit vector in the direction of p and Si.2(p) are the scalar functions, which 
depend, apart from p, on e, Ri and i?2- 
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Then, the iterative procedure is as follows. At first, the external field polarizes the 
fuUerene, giving rise to an additional field in the interior. These two fields combined 
give the field inside the fuUerene Ei = (1 — 2;)Eo, which induces a dipole moment in the 
atom Di = Q;(ti;)Ei. In the next iterative step, the dipole moment of the atom acts on 
the fullerene (problem 2), bringing about another field E2 = s(p) Di at the atom. An 
additional dipole moment is induced in the atom: D2 = a{uj) E2. For the i*^ iteration: 
Ej = s{p) Dj_i and Dj = a{uj) Ej. The total field is therefore: 

Etot(/o) = El + E2 + E3 + . . . (9) 

and the total dipole moment of the atom is Dtot(p) = «(^i^)Etot(p)- 

Evaluating the series from ([9]), one derives the following expression for the total 
field at the position of the atom: 



Etot(p,c^) ^ ^ 



Eo /^"^"j^'^'w.J nEoln 



1 + si a{uj)/Rl n 1 + (si + S2) a{uj)/Rl ' ' " ' 
The quantity z{uj) characterizes the dynamical response of the fullerene to the external 
field and is defined as: 

^ 2(iV/2 - Nn) / 1 1 \ 

i?2 ycj^ — cj^^ + iViiLU uj'^ — ujIi + 1^21^} J 1=1 
The surface plasmon frequencies are denoted by uu and UI21 for a given multipole 
moment /. These pairs of eigenmodes arise from the fullerene having a finite thickness 
and therefore two surface charge densities. The eigenfrequency ujn characterises the 
symmetric mode, in which the two surface charge densities oscillate in phase, whereas 
UJ21 characterises the anti-symmetric mode, in which the oscillations are in anti-phase [1]. 
These quantities are defined as: 

where 

pi = v/r+4i(7Ti)e^, (13) 

with ^ = R1/R2 < 1 being the ratio of the fuUerene's inner and outer radii. 
The corresponding widths and oscillator strengths of these modes are denoted by 
{Tji, Nji)j=i^2- Within the framework of our approach, the widths of these surface 
plasmon modes are not computed but are treated as parameters. 

In (ITOl) . both Si and S2 are related to the multipole expansion of the electric field in 
the situation, where a point-like dipole is positioned at p inside the fullerene: si is due 
to the field component parallel to the point like-dipole while S2 is due to the component 
perpendicular to the dipole. Their explicit forms are: 

AH ^3^. (14) 
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The quantity /3;(u;) characterizes the dynamical response of the shell to an electric 
field of the multipolarity / and the frequency uj: 

Ni2 . Nn 



UJ 



2 - + iTuLU 



+ 



UJ. 



21 



(15) 



Combining equations ([7]) and (ITU]) , one arrives at the following expression for the 
dynamical screening factor: 

1 - z{uo) ^ 

(16) 



si«(a;)//??r 



X 



S2a{uj)/Rl 



1 + (si + S2)a{uj)/Iif 



2Re 



S2a{uj)/Iif 



1 + (si + S2)a{uj)/Rl 



cos^e 



where ^ is the angle between p and Eq (see figure [T]). 



Let us apply (fT6|) to the special case where the atom is located at the central 
position. Because of the symmetry of this system, the multipole expansion of the electric 
field inside the fuUerene is limited to the dipole term. Thus, the terms Si and S2 reduce 
to: 



s,{0,u) = 2P,{lj), S2(0,^) 

and the total field at the atom simplifies to 

1 - z{uj) 



Etot(0,^) 



0, 



En. 



:i7) 



l + 2(3i{uj)a{uj)/Rl "° 
The dynamical screening factor becomes 

m-) = — (19) 

In the limiting case, where the interaction between the polarized atom and the polarized 
fuUerene shell is negligible, the denominator of equation (fT9|) goes to 1, and one obtains 



J^(0,w) 



11 - z(uj) 



(20) 



This result for JF initially derived in [1] , is independent of the nature of the endohedral 
atom. 



3. Numerical Results 

The dynamical screening factor is sensitive to a number of quantities such as the size 
and thickness of the fullerene cage, the type of endohedral atom and its position within 
the fullerene. To illustrate these dependencies, we begin by considering the case where 
the atom is located at the centre of the fullerene. Whenever possible, we compare our 
results with those presented in [2,3], where the problem of dynamical screening was 
considered within the framework of the time-dependent LDA (TDLDA). 
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3.1. Centrally positioned atom 

To investigate the impact of the type of endohedral atom on the dynamical screening 
factor, two endohedral systems are considered: Ar^lCeo and MgQCeo- In this section 
the analysis is carried out for the central position of the atom, i.e. p = (see figure [1]). 

The dynamical screening factors (fT9l) for these two systems are shown in figure [2] as 
smooth solid lines of medium thickness. For each system, it begins above the ionization 
threshold, J, of the endohedral atom, which is 16 eV for argon and 8 eV for magnesium. 
For uj < I the calculation of the screening factor according to the definition ([1]) is 
meaningless, since the atomic photoionization cross section is identically equal to zero. 

For each atom, the curve is compared with the limiting case where the interaction 
between the atom and fuUerene is neglected fl20|) . The corresponding universal curve, 
independent of the type of confined atom, is presented in both graphs as a smooth dashed 
line of medium thickness. It is seen that the screening factor can depend strongly on 
the type of endohedral atom. 

For argon, there are noticeable differences between the two cases, whereas for 
magnesium the two curves lie almost on top of each other. This can be understood 
by considering the magnitudes of the dynamic dipole polarizabilities of the atoms, given 
in figure O The polarizability of argon is large in this energy range, so the inverse effect 
of the polarized atom on the fuUerene is significant. 

In these calculations, the following parameters were used: the two dipole surface 
plasmon resonances of the fuUerene were set at 16.5 and 38.0 eV respectively, with their 
widths being 3.5 and 9.0 eV. Additionally, the mean radius and the thickness of the 
fuUerene were set to 3.54 A and 1.5 A, respectively. These parameters are the same as 
those used by Chakraborty and co-workers in the TDLDA calculations of the dynamical 
screening factor for ArQCgo [2] and MgQCgo [3]. Their results are shown in figure [2] as 
thin spiky solid lines. 

In the case of ArQCgo, there is a good agreement between the model results and the 
TDLDA calculation. The inclusion of the interaction between the polarized atom and 
fuUerene into the dynamical screening factor improves the correspondence. However, 
in the case of magnesium, there is only a qualitative agreement in the higher energy 
range. At about 10 eV, there is a strong narrow peak in the TDLDA calculation that 
cannot be described by our model. There is also a large difference about the symmetric 
plasmon mode at 16.5 eV. 

The following is a possible explanation for the large discrepancies for MgQiCeo- 
There are two types of valence electrons of the fuUerene: the a and the tt electrons [28]. 
The TT electron is more strongly bound to the carbon core whereas the a electron is more 
weakly bound. From this arises two types of surface plasmons: the tt plasmon and the 
a plasmon (see, e.g. [21,29]). The a plasmon involves mainly the more loosely bound 
electrons and is located at about 20 eV while the vr plasmon involves only the strongly 
bound ones, resulting in a lower plasmon energy of around 10 eV [21,23,30]. 



Dynamical screening of an endohedral atom 



8 




10 20 30 40 50 60 



photon energy, co, eV 

Figure 2. The dynamical screening factor for the centrally positioned endohedral 
atom: ArQCgo (top) and MgCSCeo (bottom). The common legend is given in the top 
panel. The dynamical screening factor was calculated for two cases: from (|20p . which 
neglects the atomic feedback, labelled as 'no a{u!)\ and from (|19p . with accounting 
for this feedback, labelled as 'with a(u))\ represented respectively by the solid and the 
dashed lines. These are compared with TDLDA calculations (thin solid line) from [2] 
and [3] . The ionization potential of the endohedral atom is indicated by the dotted line. 
The thick solid and dashed lines represent the factor calculated within the modified 
model, which includes the contributions of the a and it plasmons, see (f22|) and (|2T|) . The 
symmetric (s) and antisymmetric (a) modes of the tt and the a plasmons are indicated 
by arrows on the plot for MgQCeo- Regions of particularly strong interference ('int.'), 
all negative, between the tt and the a plasmons are also indicated. 
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Figure 3. Left: the dynamic dipole polarizability of argon, calculated within the 
random phase approximation with exchange (RPAE). The common legend is given 
in the right panel. Right: the dynamic dipole polarizability of magnesium, based on 
the scattering factors of [27]. The real was calculated from the imaginary part of the 
polarizability via the dispersion relation. In both plots, the solid blue line is the real 
part and the dashed green line is the imaginary part. The first ionization threshold of 
the atom is indicated by the dotted line. 



Let US introduce the two types of plasmons into the dynamical screening factor. 
To do so, one models the fullerene as a system of two co-centric spherical shells. The 
valence electrons of the fullerene are distributed into these shells, so that one shell 
contains those, which participate in the vr plasmon, and the other contains the electrons 
from the a plasmon. In reality, the plasmons will influence each other's resonance energy. 
We take the simplest approach, in which this effect is ignored. 

The finite thickness of the shells leads to the splitting of each plasmon into 
a symmetric and an antisymmetric mode. As observed in [31, 32] and seen in 
other theoretical works on these plasmon modes of the fullerene (e.g. [22, 33]), the 
antisymmetric mode does not generally manifest itself as strongly as the other mode. 

The dynamical screening factor of such a system is calculated using the iterative 
scheme described in section [2l The only modification to the scheme is the necessity of 
accounting for the two shells, which contribute to the total electric field at the centre. 
Thus, at each iterative step, there is a contribution to the electric field at the atom due 
to the vr shell and to the a shell. 

With this modification, the dynamical screening for the centrally positioned atom 
becomes (cf. (IT9l)): 

2 

1 - ^^(a;) - Za{uj) 



(21) 



The properties of the endohedral atom are encoded in its dynamic dipole polarizability 
a{u). For the case where the polarizability of the confined atom is negligible, this is 
simply 



J^iuj) ^ \1 - z^{iu) - z^{uj)\\ (22) 
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Table 1. Given below are the parameters used to calculate the modified dynamical 
screening factor (PT|) for ArQCgo and Mg@C6o- The fuUerene is modelled as two 
spherical shells - the tt shell and the a shell. Their radii are given in Angstroms. Also 
given (in cV), are the resulting plasmon energies and the widths that are assigned to 
the plasmons. 





N 


Ri 


i?2 


UJl 


UJ2 






TT shell 


80 


2.25 


3.75 


12.0 


22.2 


1 


4 


a shell 


160 


2.7 


3.7 


16.6 


37.7 


3.5 


9 



The modified dynamical screening factor is sliown in figure [2] as thick hnes. The 
sohd hne is the dynamical screening factor for the limiting case and the dashed line 
accounts for the interaction between the polarized atom and the polarized fullerene. 

For Mg@C6o, the pair of curves lie on top of each other, owing to the small dynamic 
dipole polarizability of magnesium in this energy range (see figure [3l right panel). There 
is an improvement in the correspondence between our model (12T]) and the TDLDA 
calculations of [3]. We have now reproduced the sharp peak near 10 eV that was missing 
from the original model. Within the modified model, this feature is due to the symmetric 
mode of the tt plasmon. Also reproduced is the additional feature near 20 eV arising from 
the antisymmetric mode. The positions of the 4 plasmon modes of the modified model 
are all indicated in figure [2] in the bottom panel. As most clearly demonstrated in ( l22l) . 
the screening factor consists of direct squares and cross multiplication terms, which 
raise or lower it depending on the photon energy. These cross terms, or interference 
terms, are particularly significant in 3 regions indicated in the lower panel of figure [2] by 
arrows labelled 'int.'. In these regions, the interference is strongly negative, resulting in 
well defined minima in the profile of the screening factor. There is still a quantitative 
difference near the symmetric mode of the a plasmon. It is possible to reduce this 
difference by increasing N^, shifting u;^^2 towards uJa,ii and reducing the corresponding 
widths. 

For ArQCeo, the modification of the screening factor due to the vr and the a 
plasmons has resulted in an even better correspondence with the TDLDA calculation 
from [2]. 

In these calculations, N.,^ and N^j were fixed to 80 and 160 respectively. The inner 
and outer radii and thicknesses of the two shells were adjusted to fit the energies of the 
plasmons observed in the TDLDA calculation for MgQCgo- The parameters are given in 
tabled] along with the energies uj\^2 and widths Fi 2 of the symmetric and antisymmetric 
plasmon modes of each shell. 

3.2. Arbitrarily positioned atom 

In reality, the endohedral atom is not fixed at the centre of the fullerene. Thermal 
vibrations, van der Waals interaction with the fullerene, and other effects, such as 
electron transfer, can cause the atom to be preferentially displaced from the centre. 
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Table 2. A table of the parameters of Cgo and C240 used for evaluating (fT6|l . The 
various radii of the fuUerenes are given in Angstroms, with the thickness of each 
fuUerene set to 1.5 A. The plasmon energies (in eV) are for the dipole surface plasmon 
mode only. 





R 


Ri 


i?2 


i - R1/R2 




UJ2 


Ceo 
C240 


3.5 [35] 
7.1 [34] 


2.75 
6.35 


4.25 
7.85 


0.65 
0.81 


16.9 
12.8 


33.5 
35.0 



We illustrate the spatial dependence of the screening factor with two endohedral 
species ArtQlCeo and Ar@C24o- As discussed in the previous section, there is a good 
quantitative agreement for the central case for ArCeo between the unmodified model 
f fT6|) and the TDLDA calculation of [2]. Therefore, there is no need to distinguish 
between the two types of delocalised electrons and the fullerene is modelled as a single 
spherical shell. The factor for an arbitrarily positioned atom is calculated with (fT6l) for 
argon encapsulated inside Ceo and C240 and is presented in figures H] and O respectively. 

In each figure there are 4 panels showing the screening factor (1161) as a function 
of photon energy, cu, and radial distance from the centre p (see figure [1]). Each graph 
corresponds to the indicated value of the angle 9 between p and Eq. 

The mean radius R of the C240 cage is 7.1 A [34], just over twice that of the Cgo 
cage (3.5 A [35]). The thickness of both fullerenes were set to 1.5 A. The plasmon 
energies are calculated from (IT^ . In tabled the dipole surface plasmon energies (/ = 1) 
are given together with other relevant quantities that are used in this investigation. The 
widths of all the plasmons are parameterized as Tji = 0.25ujji,j = 1, 2. 

The size of the fullerene cage has a significant impact on the dynamical screening 
factor. Firstly, as discussed in [1], the size of the cage determines the positions of the 
two plasmon modes. With a larger shell, the two plasmon modes are positioned further 
apart, which in turn controls the corresponding widths. This leads to differences in the 
energy dependence of the screening factors. Secondly, the radius of the fullerene shell, 
together with the dynamic dipole polarizability of the atom, determines the strength 
of the interaction between the atom and fullerene and how much the screening factor 
is altered. For ArOCgo, this effect is strong because of the proximity of the atom to 
the fullerene cage. For Ar@C24o, the much larger radius renders this effect negligible. 
It is this second case that causes the most prominent differences between the screening 
factors of the two systems. 

From the sets of graphs it can be seen that there is a rather weak dependence on 
6. The radial distance from the centre, p, has a much larger influence on the dynamical 
screening factor, particularly at large distances. As the atom is moved away from 
the centre of the fullerene, the screening factor (in the vicinity of the dipole symmetric 
plasmon mode), increases. This increase becomes more dramatic as the atom approaches 
the shell, as can be seen in figure |H In figure O p is only shown up to 4.3 A. In this 
region, one can see the beginning of a steep rise in the magnitude of the screening 
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Figure 4. The dynamical screening factor for ArQCeo as a function of plioton energy 
Lo and the radial distance of the atom from the fuUerene's centre, p, calculated for 
different values of the angle 6 between p and Eq. 



factor. Near the antisymmetric plasmon mode, the screening factor also increases, but 
only slightly. In the energy range between these two modes, there are small decreases 
as the atom is moved away from the centre. On the other hand, near the central region, 
the screening factor does not vary greatly. 

3.3. Spatial averaging for the screening factor 

As already mentioned, the optimal position of the endohedral atom is not necessarily 
at the centre of the fullerene. Additionally, due to thermal motion the atom is not 
fixed at any one position inside the fullerene. For meaningful observable results for the 
dynamical screening factor, it is necessary to perform an averaging over all possible 
positions of the endohedral atom. 

A spatial average of the dynamical screening factor {T{uj)) can be carried out as 
follows: 

(^H)= f F{p,uj)W^{p)dp. (23) 
Jv 

This integration is carried out over the volume V where the atom may reside: a sphere 
of radius Ri — Ra, where Ra is related to the size of the confined atom. When the atom 
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Figure 5. Same as in figure H] but for Ar@C24o- 
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is sufficiently close to the fuUerene cage, quantum effects such as hybridization become 
significant and the suitability of this semi-classical model is reduced. 

The Boltzmann distribution W^{p), governing the probability of finding the atom 
at p and temperature T, is given by: 

f/tot(p)-f/tot(0)' 



Wt(p) = A exp 



kT 



(24) 



where k is the Boltzmann factor and A is the normalisation constant obtained from 

The total interaction between the fuUerene and the confined atom can be presented 
as a sum of two terms: 

f/tot(p) = t/(p) + f/vdw(p). (25) 
The term, U{p), is the potential energy of the interaction of the atomic dipole moment 
D(p) with the total electric field Etot(p): 

1, 



Uip) 



-Re(E:„,(p)D(p)) 



(26) 



The second term in (p5l) . f/vdw? stands for the van der Waals interaction between 
the fuUerene and the endohedral atom. This is approximated by the Lennard- Jones 
potential, and is given by: 

Cl2 Cq 



\12 



(27) 
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Figure 6. The f/vdw profile for ArQCeo (left) and Ar@C24o (right) are shown. 



where the sum is carried out over all atoms of the fuUerene and rj is position of the 
i-th atom. The parameters for the interaction between an atom of the fullerene and 
the argon atom are taken from [36]: the depth of the potential well is 0.005 eV and the 
equilibrium distance between the two atoms is 3.84 A. This translates to Cg = 53.8 a.u. 
and Ci2 = 3.93 x 10^ a.u. respectively. 

For the evaluation of the van der Waals interaction energy, the fullerene is 
approximated as an infinitely thin sphere of radius R and the sum of equation ( 1271) 
becomes: 

f/vdw(p) = / - 5) (28) 

s 

The integral is carried out over the surface of the sphere, eliminating the angular 
dependence of van der Waals interaction energy. The distance between the surface 
element dS and the endohedral atom is denoted by a. The quantities and ai2 for the 
fullerene of n carbon atoms are defined as 

nCa nCi2 , . 

= 4^' = 4^' ^2^) 

where u/AttE? is the surface density of carbon atoms. Evaluating the integrals, one finds 
that the van der Waals interaction energy between the fullerene and the endohedral atom 
is 



where 



^vdw(p)=t/l2(p)-f/6(p), (30) 

2TxRaQ f I 1 




4p \{R~pY {R + p)\ 
2TiRai2 ( \ 1 ^ ^ ' 



lOp V(^-p)'° (^ + P)'°, 

The van der Waals potential profile for the two endohedral systems are shown in figure 

El 
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In both cases the dominant term in the total potential energy is the van der Waals 
term U^dwip)- ^^e case of Ar@C6o, the radius of the fullerene (3.5 A) is less than 
the equilibrium distance of the C— Ar van der Waals interaction (3.84 A). With the 
exception of a small central region (p < 0.6 A) the repulsive term of the Lennard- Jones 
potential dominates, as shown in the left panel of figure El The atom is, therefore, 
strongly confined to the centre of the fullerene. For Ar@C24o, the size of the cage (7.1 
A) is sufficiently large that there is a potential minimum, which lies in a spherical shell 
of radius of 3.6 A. The depth of this potential well is 0.1 eV (see right panel of figure 

ED- 

The results of the spatial averaging are presented in figure [71 In all cases, the 
external field strength was set to 0.001 a.u. (where 1 a.u. of electric field strength is 
5 X 10^ V/cm). For ArtQiCeo, the argon atom is essentially fixed at the centre. The 
results of the spatial averaging refiects this. Large increases in temperature only allows 
a minimal increase in the mobility of the endohedral atom, and has almost no impact 
on the spatially averaged dynamical screening factor. For Ar@C24o, the argon atom 
should have enough thermal energy at T ^ lOOO/C to escape the potential well and 
reach the central region of the fullerene. However, no significant changes are observed 
in the spatially averaged dynamical screening factor as the temperature is changed from 
lOOK up to 1500ir. 

4. Summary 

In this semi-classical approach, a dielectric formalism was taken to study the dynamical 
screening of the endohedral atom. The mutual interaction between the polarized atom 
and the polarized fullerene is now accounted for. Previous work [1] has been generalised 
to consider an atom confined at an arbitrary position inside the fullerene. 

A more indepth study of the central case was presented, showing the effect of 
accounting for the polarizability of the endohedral atom. Case studies were performed 
on ArQCgo and MgQCeo- We show that the dynamical screening factor can depend 
strongly on the type of endohedral atom. 

A comparison of the dynamical screening factor for ArtQCgo and MgtOlCeo (the atom 
was positioned at the central position) was performed with the TDLDA calculations 
[2,3]. A good agreement was found for argon, however, only a qualitative agreement 
was observed for magnesium. 

In an attempt to explain the discrepancies between our model and the TDLDA 
calculation, the dynamical screening factor was modified to include the two types of 
plasmons of the fullerene, the n and the a plasmons. This was achieved by modelling 
the fullerene as a system of two co-centric spherical shells, each being associated with 
one of the types of plasmons. The modified model shows an improved correspondence 
with the TDLDA results. 

Through the endohedral species ArCiCgo and Ar@C24o, we show the spatial 
dependence of the dynamical screening factor. There is only a weak dependence on the 
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Figure 7. The spatially averaged dynamical screening factor for ArQCeo (top) and 
for Ar@C24o (bottom) for T lOOA' (thick sohd hnc) and T = 1500X (think dashed 
line). The common legend is shown in the top panel. The screening factor for the 
centrally positioned atom is shown for comparison (thin solid line). 



angle between p and Eq. The distance between the endohedral atom and the fuUerene 
cage has a significant impact on the screening factor - as the atom approaches the cage, 
its magnitude increases rapidly. 

The size of the cage also impacts the energy dependence of the screening factor. The 
response of the fuUerene to the external field - the resonance energies of the symmetric 
and antisymmetric dipole plasmon modes - is controlled in part by its radius. This was 
discussed in [1]. In the present paper, the radius of the shell has an additional effect - it 
determines how strongly the screening factor is modified by the interaction between the 
atom and the fuUerene. For larger fullerenes, this interaction is diluted by the increased 
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distance between the atom and the cage. 

As the atom is not fixed to any position within the fullerene, a Boltzmann 
distribution was chosen to govern the probabihty of the confined atom being located at 
some arbitrary position within the fullerene. This distribution is based on the interaction 
potential energy (electrostatic and van der Waals) between endohedral and the cage. 
For the two cases studied, the van der Waals interaction is the dominant term. The 
argon atom is of similar size to the Ceo cage and is therefore located at the centre. In 
the case of C240, the cage is sufficiently large for the atom to move freely within it. The 
van der Waals interaction between argon and the cage confines the atom to a spherical 
shell of radius of 3.6 A. In both cases, the spatially averaged dynamical screening factor 
is only weakly dependent on the temperature. 

Despite the simplicity of this model, it provides an instructive description of the 
phenomena discussed here. This model may be further applied to study the dynamic 
responses of more complicated situations such as multi-shelled fuUerenes. 
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